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DIGEST 


The  principles  governing  freeze-drying  rates  of  heat  and  mass  trzmsfer 
have  been  reviewed  and  defined  in  practical  terms. 

The  various  methods  and  resistances  to  heat  treinsfer  are  discussed 
eind  practical  use  is  made  of  the  equations  presented. 

Mass  transfer  through  the  various  phase  changes  is  discussed.  Theo¬ 
retical  discussions  of  the  mode  of  mass  transfer  at  each  step  and  suitable 
relationships  are  developed  which  describe  the  limiting  conditions  for  the 
mass  transfer  rates-  By  use  of  the  over-all  mass  transfer  equations,  rates 
of  drying  and  total  times  required  may  be  determined. 

Sample  calculations  illustrating  the  use  of  the  relationships  and 
theoretical  considerations  are  included. 
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I.  INTRODUCTION 


Freeze-diying  is  a  method  of  drying  materials  by  sublimation,  using 
low  temperatwes  and  reduced  pressures.  It  is  generally  employed  where 
heat-sensitive  materials  must  be  processed  to  a  dry  state.  The  principle 
used  is  one  in  which  the  material  to  be  dried  is  frozen  in  some  convenient 
form;  slabs  or  pellets  are  most  frequently  used.  The  frozen  material  is 
then  subjected  to  a  high  vacuum  which  is  below  the  vapor  pressure  of  water 
at  the  temperature  of  the  material  to  be  sublimed.  Under  these  conditions, 
the  ice  (or  frozen  material)  sublimes  from  the  frozen  mass  amd  is  collected 
on  a  cooled  condenser  whose  temperature  is  such  that  the  vapor  is  desub- 
limed.  As  sublimation  continues,  heat  is  withdrawn  from  the  material  being 
dried.  In  order  to  maintain  the  temperature  of  the  material  high  enough  so 
that  the  sublimation  process  will  not  stop,  heat  is  supplied  to  the  material 
from  an  evaporator.  The  rate  at  which  heat  is  supplied  is  dependent  on  the 
drying  cycle  desired. 

A  typical  apparatus  consists  of  a  horizontal,  cylindrical  body  contain¬ 
ing  a  horizontal  evaporator  zmd  condenser.  The  evaporator  and  condenser 
consist  of  plates  or  coils  in  which  a  suitable  heating  and  cooling  medium 
may  be  circulated.  Cast  metal  with  a  high  thermal  conductivity  and  a 
blackened  surface  to  increase  transfer  of  heat  by  radiation  is  used  as 
materizd  of  construction  for  the  trays  which  hold  the  product  to  be  dried. 
The  material  trays  may  have  ribs  cast  as  an  integral  part,  to  provide  more 
uniform  heat  transfer  through  the  layer  of  material  being  dried.  Sources 
of  refrigeration,  heat,  and  vacuum  are  connected  to  the  dryer  in  order  to 
pro\’ide  proper  progrzimming  of  heat  transfer. 

Figure  1  shows  a  typlw^l  in^-tallation,  including  the  arrangemen-^  of  the 
basic  components  of  the  batch  system.  This  report  covers  the  theory  of 
freeze-drying  by  the  batch  method,  but  does  not  cover  tteory  specific  to 
continuous  systems  such  as  that  described  by  Maister.-/* 


See  Literature  Cited. 
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Fi^re  1.  Equipment  Arrar  geraent  in  Typical  Batch  Freeze-Drying  System 
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II.  THEORY  OF  HEAT  AND  HASS  TRANSFER 


A.  GENERAL  CONSIDERATIONS 

The  process  of  heat  transfer  in  freeze-drying  can  be  stated  as  follows. 
Heat  must  be  transferred  from  the  heater,  throu^  the  evaporator  shell, 
through  the  vapor  space,  through  the  frozen  material  to  the  interface 
where  sublimation  is  actually  taking  place.  Heat  is  also  transferred  from 
the  desubliming  vapor  to  the  ice  surface,  through  the  ice  cake  and  the  con¬ 
denser  wall  to  the  refrigerant,  idiere  it  is  absorbed  and  carried  off. 

The  process  of  mass  tramsfer  begins  at  the  subliming  interface.  The 
vapor  formed  there  passes  tlirough  the  porous  cake  of  previously  dried 
material  and  into  the  vapor  space  separating  the  material  and  the  condenser o 
From  here  it  passes  to  the  interface  between  the  vapor  space  and  the  ice 
cake  where  it  is  desublimed  and  the  process  of  mass  transfer  is  ended. 

It  is  these  several  steps  of  heat  and  mass  transfer  that  actually 
govern  the  rates  of  drying, 

1,  Drjdng  Stages 


There  are  two  stages  in  the  drying  process:  the  constant-rate 
period  and  the  falling-rate  period.  The  point  >rtiich  separates  these  two 
periods  is  termed  the  critical  moisture  content  (CMC),  Experimentally, 
this  point  is  determined  in  the  drying  process  when  the  temperature  of 
the  drying  material  begins  to  rise,  because  the  constant-rate  period  is 
characterized  by  a  constant  temperature  of  the  material.  The  CMC  repre¬ 
sents  the  point  where  all  of  the  free  moisture  has  been  removed  from  the 
material. 


Molochnaya^  has  proposed  the  following  equation  for  the  moisture 
content  at  this  point: 


^  0.5  ^  0.75 

A  (tm)  (h) 


(1) 


where  A  =  empirical  constant  (for  gelatin,  6410) 
2h  =  thickness  of  specimen  being  dried 
Pm  =  pressure  in  system 
tm  =  final  temperature  of  specimen 
w^  =  CMC 

w 


(2) 


Laabert^/  gives  this  equation; 
w  = 

"  As 

where  As  =  heat  of  sublimation 

Eq  =  (353) -(lb  water  sublimed) /(total  watt-hr  used) 

During  the  constemt-rate  period,  heat  is  supplied  at  a  constant 
rate  to  the  dryer.  The  faster  heat  can  be  supplied  under  the  limitations 
of  the  heat  and  mass  transfer  rates,  the  faster  the  drying  process.  This 
heat  is  transferred  thi*ough  the  cake  and  used  to  sublime  the  ice.  All  of 
this  added  heat  is  used  in  sublimation,  allowing  the  material  to  remain 
at  a  constant  temperature.  Once  the  CMC  is  reached,  the  temperature  of 
the  material  begins  to  rise  because  all  of  the  free  moisture  has  sublimed, 
resulting  in  less  evaporative  cooling.  At  this  time,  the  heat  to  the  dryer 
must  be  reduced  to  prevent  overheating  of  the  material  and  possible  melting 
of  some  portions.  The  amovint  of  heat  reduction  is  dependent  on  the  rates 
of  heat  and  mass  transfer.  At  the  firjt  evidence  of  temperature  increase 
at  any  point  of  the  material  the  heat  should  be  reduced  to  avoid  toasting 
at  anj'-  one  point,  even  though  there  may  be  free  water  in  other 

portions  of  the  material. 

Methods  of  heating  include  resistance  wires,  induction  heating, 
high-frequency  fields,  dielectric  heating,  and  infrared  lamps. 

In  I'losdorf’s  bookl/  Greaves  points  out  that  the  control  of  the 
supply  of  latent  heat  of  sublimation  is  the  most  important  single  factor  in 
freeze-drying  rates.  How  quickly  heat  can  be  applied  will  depend  upon: 

(a)  Capacity  of  refrigeration  at  condenser  temperature  reading; 

(b)  Highest  safe  temperature  at  which  material  may  be  dried  (for 

a  given  heat  input,  thi*}  temperature  depends  upon  the  degree  of  obstruction 
to  vapor  flow  and  upon  condenser  temperature) ; 

# 

(c)  Highest  temperature  at  which  it  is  safe  to  operate  heaters;  and 

(d)  9.ate  of  transfer  of  heat  through  frozen  material. 

The  heat  applied  to  the  frozen  material  during  drying  produces  a  vapor 
pressure  difference  between  the  evaporating  and  condensing  sixrfaces.  This 
reflects  a  difference  in  temperature  between  the  two  surfaces. 

Ssjice  rate  of  flow  depend?  upon  heat  input  (ClapejTon's  equation), 
a  fommla  may  be  written: 


(-3) 


Rate  of  flow  = 


_ Ap _ 

obstructive  resistance 


-  (K)  * 


(Fatts) 


where  Ap  ^  vapor  pressure  difference 


K  =  voltme  constant  in  watt-hours  per  milliliter 

Greaves£/  defines  the  unit  of  obstructive  resistance  siS  ttiat  which, 
under  a  vapor  pressure  difference  of  0^01  mm  Hg,  passes  vapor  at  the  rate 
at  idiich  it  is  liberated  by  a  heat  input  of  one  watt.  Therefore; 


I  _JL 
R  “  ioo 


(4) 


wiere  P  =  vapor  pressure  across  the  resistance  in  mm  Hg 
R  =  obstructive  resistance 
ff  =  heat  input  in  watts 
2 .  Optimum  Cycle 

During  the  first  stage  of  drying,  heat  must  be  introduced  into  the 
frozen  material  as  rapidly  as  possible  without  causing  it  to  soften  or 
melt.  At  the  same  time,  a  maximum  rate  of  flow  of  water  vapor  away  from 
the  evaporating  surface  mtist  be  established.  To  accomplish  this  rapid 
vapor  flow,  adequate  passageways  must  be  provided  for  the  vapor,  which 
must  then  be  condensed  or  evacuated  efficiently.  However,  rupply  of ^eat 
to  frozen  material  is  the  controlling  factor,  according  to  Greaves.^ 

The  basic  limitaiiosis  listed  by  Greaves^./  include; 

(a)  Heat  must  not  be  carried  to  the  walls  of  the  container  >diich 
holds  the  product  faster  than  the  heat  can  be  conducted  through  the  frozen 
mass  to  a  free  surface  where  it  is  used  to  induce  sublimation.  Otherwise, 
melting  id.ll  occtir  adjacent  to  the  container  wall. 

(b)  Heat  can  be  carried  down  directly  to  the  evaporating  surface 
to  avoid  conducteince  through  ice,  but  even  here  there  is  a  limitation.  As 
soon  as  sublimation  has  proceeded  to  an  appreciable  extent,  the  ice  layer 
has  receded  from  the  level  of  the  original  surface  so  that  the  evaporating 
surface  is  confined  within  the  interstices  of  the  outer  frame-work  of 
porous  dry  solid.  The  heat  must  then  be  carried  across  this  porous  struc¬ 
ture,  but  the  temperature  of  this  dry  portion  of  the  material  must  not  be 
brought  up  to  the  level  at  which  it  is  harmed. 


Flosdorfl/  has  detenained  that  the  difference  between  the  partial 
pressure  of  water  vapor  that  can  be  established  within  the  vacuum  system 
the  condenser  and  the  vapor  pressure  of  the  evaporating  material  determines 
the  theoretical  maximum  flow  rate  of  vapor.  The  vapor  pressure  of  the 
material  at  the  temperature  at  trfiich  it  should  be  dried  determines  Ihe 
temperature  at  which  the  condenser  should  be  maintained.  Differential 
pressure  between  the  partial  pressure  and  the  vapor  pressure  is  the 
driving  force. 

For  each  of  several  temperatures  there  is  a  m-in-imai  temperature  of 
condensing  surface,  below  which  there  is  very  little  further  increase  by 
further  reduction  in  condenser  temperature.  To  illustrate  this  point,  x 
on  each  ctirve  in  Figure  2  represents  the  condenser  (ice)  surface  tempera¬ 
ture  giving  rise  to  maximum  differential  according  to  the  Napier  equation 
(laws  of  adiabatic  gas  flow  through  orifices),  assuming  no  restrictive 
orifices  in  pipe  lines,  etc.,  other  than  in  the  interstices  of  the  porous, 
dry,  outer  layer  of  material  itself.  This  applies  to  the  first  stage  of 
dr^ng  when  ice  is  being  sublimed. 

Farther,  as  drying  proceeds,  it  is  limited  in  rate  by  diffusion  of 
vapor  through  the  interstices  of  the  porous  dry  outer  layer  of  material, 
which  act  as  orifices.  By  the  Napier  equation  a  differential  between  the 
^'apor  pressure  at  the  surface  of  the  condenser  and  the  vapor  presstare  at 
the  ice  surface  within  the  material,  where  the  former  is  55  per  cent  of  the 
latter,  will  result  in  the  maximum  rate  of  flow  obtainable.  For  example, 
this  condition  is  satisfied  by  a  material  temperature  of  -18°C  (for  ice) 
and  a  condenser  temperature  of  -25°C.  The  drying  time  is  21  hoiu^. 

During  another  run  with  a  condenser  temperature  of  -40°C,  the  dry¬ 
ing  time  was  the  same  but  the  jacket  temperature  of  the  drying  chamber  had 
to  be  raised  from  60  C  to  9C°C  in  order  to  keep  the  material  at  -i8°C.  Tf 
the  jacket  temperature  had  not  been  raised,  drying  would  have  been  retarded 
to  27  hours,  with  a  -40°C  condenser  and  the  material  at  ..bout  -32°C.  The 
presstire  was  500  microns  in  both  caser;,  but  if  the  jacket  had  not  been 
raised  to  90°C  with  a  -40®C  condenser,  this  pressure  would  have  had  to  be 
decreased  to  225  microns.  The  efficiency  of  transfer  of  heat  to  the 
material  is  all-important  in  rapid  drying,  provided  that  the  ice-condensing 
surface  is  below  the  critical  level  discussed  above  for  obtaining  the  most 
rapid  flow  of  vapor. 

Reduction  of  the  vapor  pressure  on  the  condensing  stirface  to  less 
thzin  55  per  cent  of  the  pressure  on  the  drying  side  of  each  orifice,  by 
laws  of  adiabatic  gaseous  flow,  can  produce  no  further  increase  in  rate 
of  flow  tnrough  the  orifice.  In  other  words,  beyond  a  certain  point  lower¬ 
ing  the  temperature  of  the  condenser  cannot  compensate  for  small  vapor 
lines  or  restri.ctions . 


Pressure  Difference,  microns  of  mercury 


Figure  2.  Driving  Force  in  Terms  of  Differential  in  Pressure  Produced  by 
Various  Condenser  Temperatures  for  Products  Maintained  at  Four 
Different  Temperatures. 


By  use  of  condensers  of  great  surface  having  constant  efficiency 
and  by  a  proper  balance  of  all  other  factors  concerned,  a  very  important 
practical  conclusion  is  obvious  as  a  result  of  applying  the  law  of  flow» 

It  becomes  necessary  to  accept  an  ice  surface  with  its  poor  heat  transfer* 
However,  by  use  of  condensers  of  large  surface  area  where  the  heat  to  be 
transferred  per  unit  area,  is  kept  small,  a  higher  surface  temperature  is 
adequate.  Thus  a  higher  temperature  of  refrigerant  can  be  permitted. 
Condensers  may  be  operated  at  temperatures  much  above  any  originally 
believed  possible  without  sacrifice  of  effectiveness.  With  material  at 
-21*^0,  a  temperature  of  ice  on  the  condenser  of  about  -27°C  produces  a 
vapor  pressure  equal  to  55  per  cent  of  that  of  the  material;  thus  a 
lower  temperature  at  the  ice  surface  is  no  advantage.  Also,  for  greatest 
efficiency,  the  temperature  of  the  material  being  dried  must  be  no  lower 
than  necessary. 

The  final  temperature  to  which  the  material  may  be  taken  in  the 
second  stage  of  drying  and  the  vapor  pressure  of  the  product  at  that 
temperature  also  have  a  bearing  on  the  required  temperature  for  the  con¬ 
denser  or  pressures  needed  in  direct  pumping.  This  is  related  to  the 
final  residual  moisture  content  obtainable.  Whatever  the  vapor  pressure 
of  the  dry  product,  the  condenser  or  pump  must  establish  a  pressure  that 
is  lower. 

In  cases  of  hygroscopic  products  not  stable  at  higher  temperatures, 
a  regenerable  desiccant,  either  calcium  sulfate  or  perchlorate,  may  be 
used  advantageously  to  yield  a  low  partial  pressure  of  water  vapor  during 
the  final  stage.  Since  such  a  small  actual  weight  of  water  is  removed  in 
this  stage,  regeneration  of  the  desiccant  is  infrequent.  This  is  a  less 
expensive  means  than  refrigeration  to  produce  lower  final  condenser 
temperature  or  a  steam  eject'^r  for  lower  total  pressure. 


B.  SUBLIMATION 

1.  General  Considerations 

J.B.  Lambesri  presents  a  very  complete  review  of  sublimation  in  his 
thesis..^/  Almost  all  of  the  material  that  follows  has  been  taken  from  this 
work.  Equations  generfxily  apnly  to  slabs,  but  with  appropriate  mathematical 
manipulation  may  be  apolied  to  any  shape. 

At  the  subliming  interface,  the  absolute  rate  of  sublimation,  w_*, 
is: 

Ws»  =CXp5*  y277'RTs* 


(5) 


n 


where  (f  =  accommodation  coefficient 

=  vapor  pressure  in  equilibrixan  with  Tg# 

Tg«  =  temperature  of  subliming  interface  in  degrees,  abs. 

MW  =  molecular  weight  of  subliming  substeuice 
R  =  gas  constant 

The  accommodation  coefficient  is  used  to  compensate  for  the  fact 
that  some  molecules  collide  with  the  interface  and  may  be  reflected  instead 
of  sublimed.  Since  the  state  of  the  gas  does  not  affect  the  total  rate  of 
vaporization,  Equation  (5)  is  used  to  calculate  evaporation  rate  under 
nonequilibrium  conditions. 

Previous  investigations  taike  the  net  sublimation  rate,  w,  equal  to 
the  difference  between  evaporation  and  condensation  rates  at  the  interface, 
calculated  by  Equation  (5)  and  assuming  temperature  equilibrium  but  not 
pressure  equilibrium  at  the  interface,  therefore: 


where  Pq  =  initial  vapor  pressure  of  gas 

The  same  accommodation  coefficient  is  used  for  evaporation  as  for  condensa¬ 
tion. 

Schrage.^^  cites  an  error  in  this  reasoning.  Equation  (5)  for  cor 
densation  assumes  that  molecular  velocity  distribution  of  gas  in  the  vapor 
space  is  Maxwellian  with  respect  to  coordinates  fixed  in  space.  Since  net 
mass  transfer  occurs,  the  system  is  not  in  equilibrium  and  ordinary  Maxwell 
distribution  does  not  apply.  Schrage  assumes  that  gas  velocity  distribu¬ 
tion  becomes  Maxwellian  with  respect  to  the  molecular  streeun  velocity  in  a 
very  few  meeui  free  paths  from  the  interface.  Hierefore: 


where  P  =  correction  factor  to  Equation  (5),  which  is  itself  a  fimction 
of  net  mass  transfer  rate, 

Tq  =  gas  temperature  at  interface,  K° 

Wq  =  absolute  condensation  rate  at  interface 
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The  factor  P is  defined  by  equations: 


0  '  27r^/2 


w  * 
"s 

I  / 


P  * 


'To\  -1/2 


"0^  -  cpir" 


/2 


(1  -  erfcy) 


fx  2 


vrhere  erfcy  =  ■  error  integral. 


(0) 

(9) 


Addition  of  total  evaporation  and  condensation  rates  from  Equations  (5) 
and  (7)  yields  net  transport: 


1/2 

=  v;-* 

3 

1- 

T 

o  , 

V  / 

r 

(10) 


The  factor  P  is  uniquely  related  to  v/ws*  for  the  special  case  >diere  Tp/Tg<^ 
is  equal  to  one.  If  eliminated  from  ^nations  (9)  and  (10), ©may 

be  plotted  as  a  function  of  w/wg*.  Values  of  p  are  then  calculated  from 
various  w/w^*  values  b^'  using  Equation  (9)  Kith  the  plot  shown  in  Figure  3. 

In  applying  Schrage'sli  theory  to  the  case  of  To/Tg*  =  1  for  net  evaporation, 
the  ratio  Pq/Ps'®’^  decreases  with  w/wg*;  it  increases  for  net  condensation. 
Schrage's  modified  theory  predicts  these  critical  or  maximum  mass  transfer 
ratios: 


-  0.  C-.-r  0Vcjp.;r't+'i  --n 

'''  ''"s""  =  13  fi-r  condensation 

The  above  theories  should  be  applied  to  mass  transfer  processes  not  far 
from  equilibrium.  Ihe  accommodation  coefficient  for  ice-v/ater  has  been 
found  to  be  one,  assvuning  no  temperature  dependency, 

.  • 

Kramers  and  StemerdingJ.-  give  equivalent  equations  for  sublimation 
rates.  On  the  basis  of  Knudsen's  formula,  they  propose  a  theoretical 
sublimation  rate,  w^: 


^'o 


M 


'2j-fR  T, 


es 


(11) 


where  M  =  molecular  weight  of  evaporating  substance 
77gg  =  vapor  pressure  of  evaporating  substemce 


Transfer  Ra»e,  (Eqoo*'®" 


t 

\ 

! 
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R  =  gas  constzmt 

Tgg  *  absolute  temperature  of  surface 

This  assumes  an  accommodation  coefficient  of  unity,  uhich  implies 
no  rebounding  of  any  of  the  molecules  on  collision. 

In  freeze-drying  practice  a  similar  formula  is  used: 


I.  M  ... 

27rR  Tes 


where  f^  the  coefficient  of  evaporation,  is  between  zero  and  one. 

In  this  formula  the  following  effects  causing  w  to  be  less  than  wq 
may  have  been  collected: 

(a)  Accommodation  coefficient  is  less  than  one. 

(b)  Pressure  in  vapor  space  is  not  zero,  so  that  vapor  molecules 
have  chance  of  rettxming  to  the  evaporator  surface.  This  pressure  will  be 
at  least  the  vapor  pressure  (TTcs)  of  the  condenser  surface,  which  is 
eventually  increased  as  a  result  of  resistance  to  vapor  flow  from  the 
evaporator  to  the  condenser  either  by  diffusion  through  noncondensable 
gas  (air)  or  by  wall  friction.  The  latter  effect  is  important  in  freeze 
drying  because  the  layer  of  dry  product  obstructs  this  flow. 

(c)  (f,  the  accomaodatio.^'  coefficient,  may  be  limited  by  the 

at  which  heat  is  being  supplied  to  the  evaporator  and  by  the  rate  of  -f.-.i'-i- 
drawal  of  the  vapor  by  condensation,  absorption,  or  pumping: 

Cr=  0.94  i  0.06  for  ice  -60°C - >  -SS^^C 

If  there  is  no  frictional  resistance  between  evaporator  and  condenser,  and 
if  the  total  water  vapor  plus  air  pressure  in  the  system  is  ^  =  (?  +ff  )» 
the  rate  of  sublimation  is  given  by: 


where  M  =  molecular  weight  of  evaporating  stirface 
R  =  gas  constant 
D  -  coefficient  of  diffusion 
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P  =  total  pressure 

Tf'  =  pai'tial  pressure  of  water  vapor 
T  =  average  absolute  temperature 

1  =  disteince  from  evaporator  surface  to  condenser  surface 

-  partial  pressure  of  evaporating  surface 

=  partial  pressure  of  condensing  surface 

The  theory  is  that  two  resistances  exist,  in  series,  to  the  transfer  of 
vapor.  The  term  Wq  (1  -  U—)  gives  the  net  transport  by  free  evaporation 

Ties 

between  the  evaporator  surface  and  the  adjacent  vapor,  which  has  partial 
water  vapor  pressuj'e  TT,’  (TTes  “  TT*  ^  7Tcs^*  right-hand  expression 
represents  the  diffusional  flow  of  water  vapor  through  air  from  Jt '  neeir 
the  evaporator  to  77cs  condenser  surface,  over  distzmce,  1,  between 

them. 


2.  Special  Cases 


1,  If  the  air  pressure  at  the  condenser  is  relatively  lor,  or 
P  =  TTruf  also7r''^P,  giving 

L5 


w  =  w. 


(14) 


2.  If  w  is  much  less  than  e.g.,  by  an  insufficient  supply  of 
heat,  it  follows  from  Equation  (13)  that  Jf'  "=  /Tes*  absence  of 

air  this  means  that  TTes  ®  TT-o  i*e.,  the  surface  temperatures  of  evaporator 
and  condenser  are  equal.  In  the  presence  of  air  the  equation  becomes: 


MDP 

Rn 


ui 


p  -  m 


cs 


P  -  TT, 


es/ 


where  the  symbols  are  the  same  as  those  in  Equation  (13). 


(15) 


3.  Application 

For  the  application  of  the  theory,  a  McLeod  vacuum  gauge  gives 
the  total  pressure,  P*H20  +  P*  air,  and  the  ice  surface  temperatures  tgg 
and  tgg  have  to  be  known.  They  can  be  derived  from  the  meastired  tempera¬ 
tures  tg  and  tg  by  correcting  for  the  temperature  drop  over  the  Iziyer  of 
air.  The  latter  quantity  can  be  calculated  from  the  heat  flow  corresponding 
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to  the  measured  rate  of  sublimation  and  the  thickness  and  thermal  con¬ 
ductivity  of  the  ice  layers.  In  experiments  carried  out  for  t^  between -2(Pc 
-2‘5°C  2ind  tp  between  -30°C  and  -37°C  at  an  air  pressure  less  than  one 
micron  of  mercury,  the  rate  was  of  the  order  lO"*^  gn/cm^sec;  i.e.,  two 
to  four  per  cent  of  the  maximum  rate  at  these  temperatures.  Consequently, 
in  these  runs  tgs  is  approximately  equal  to  t^g.  Prom  the  measured 
temperature  diff^erences ,  tg-t^,  and  the  corresponding  rates  of  evapora¬ 
tion  and  ice  thickness,  the  thermal  conductivity  was  found  to  range 
between  5,0  x  lO”^  and  6.1  x  10“^  calories  per  centimeter  degree  Centi¬ 
grade,  giving  an  average  value  of  5.7  +  0.4  x  10“^  at  a  mean  temperature 
of  -25°C,  The  agreement  with  the  generally  accepted  value  of  5,7  x  10”3 
at  -20°C  is  satisfactory. 


The  following  situation  demonstrates  the  existence  o“  a  maximum  rate 
of  evaporation.  Ice  is  evaporated  at  a  certain  constant  rate  with  a  con¬ 
denser  temperature  of  about  -30°C  ard  lowest  possible  residual  air  pressure, 
so  that  tes  =  t^g.  If  the  condenser  tempera tmre  is  lowered  and  the  sub¬ 
limation  rate  remains  the  same,  a  difference  between  tes  ‘•^cs 
appear.  It  increases  in  such  a  way  that  finally  tgs  will  become  constant 
for  very  low  condenser  temperature.  This  is  shown  in  Figure  4. 


It  is  evident  that  at  the  corresponding  temperatures  tes  '^csi 
these  are  the  maximum  rates  which  can  be  compared  with  the  equation  w  =  Wq 

— As  an  example,  the  lowest  curve  has  a  point  tgg  =  -57°C,  tgg  = 


1 


n 


es 


-64°C.  Kith  these  values,  w  gives  a  rate  of  1,3  x  10"^  gm/cm^sec,  whereas 
the  measured  rate  was  1.2  x  10"'*.  Agreement  is  good,  showing  an  accommoda¬ 
tion  coefficient  of  nearly  3.0. 


Table  I  presents  data  for  air  pressures  less  than  one  micron  of 

mercury. 


The  values  of  unity  for  Wexpt/"caic  support  the  above  theory.  This 
is  especially  important,  because  in  most  experiments  the  mean  free  path  of 
the  molecules  was  less  them  the  distances,  1,  whereas  Equation  (13)  is 
generally  accepted  only  for  cases  in  which  the  mean  free  path  greatly 
exceeds  1.  Since  Equation  (13)  has  been  derived  from  an  accommodation 
coefficient  of  unity,  these  resu3.ts  show  that  this  assumption  is  very 
probably  correct  in  the  temperature  range  from  -40  to  -60°C, 

Only  at' low  condenser  temperatures  and  at  relatively  great  distances 
between  the  evaporator  eind  condenser  are  rather  small  WgxtpA^calc  ratios 
found.  Probably,  wall  friction  and,  to  a  minor  extent,  the  pressure  of  the 
residtial  air  are  responsible. 


Temperature  of  Evaporator  Surface,  tAc  (°C) 
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TABLE  I.  SUBLBIATION  RATES  FOR  ICE  AT  PRESSURES  LESS  ITIAN 

ONE  MICRON  OF  MERCURY 


ies 

^cs»°^ 

t 

‘■es-cs*  ^ 

Calctalatedi^/ 

Sublimation 

Rate,104 

gm/cffl^sec 

Experimental 
Sublimation 
Rate.lO^ 
go/ cm^sec 

w  expt. 
w  calc. 

12 

-44,0 

-47.9 

3.9 

3.6 

5.21 

1.4  t  0.2 

-49.2 

-57.0 

7.8 

3.6 

3.54 

1.0  ±  0.15 

-52,6 

-61.8 

9.2 

2-4 

2.09 

0.9  i  0,15 

-54,8 

-65,5 

9.3 

2.0 

1.55 

0.8  ±  0.20 

-57,0 

-64.9 

7.9 

1.5 

1.16 

0.8  ±  0.15 

-58,8 

-64,4 

5.6 

0.80 

0.59 

0.7  i  0.15 

-61.7 

-70.0 

8.3 

0.75 

0.675 

0.9  i  0.05 

20 

-41,3 

-47.6 

6.3 

7.0 

6.60 

0.9  ±  0.05 

-46,2 

-52.7 

6.5 

4.4 

5.00 

1.1  i  0.1 

-46,0 

-54.6 

8.6 

5.1 

3.72 

0.7  1  0.1 

-50,4 

-58.9 

8.5 

3.1 

3.31 

1.1  i  0.1 

-50,4 

-58,6 

8.2 

3.0 

2.78 

0,9  +  0.1 

-52,2 

-62.2 

10.0 

2.6 

2.23 

0.9  t  0.1 

-54.3 

-64.2 

9.9 

2.0 

1.60 

0.8  i  0.05 

-56.4 

-66.7 

10.3 

1.6 

1.15 

0.7  1  0.05 

-61.3 

-70.1 

8-8 

1.0 

0.607 

0.6  t  0.05 

44 

-43.3 

-47.6 

4.3 

4.3 

2.70 

0.6  1  0.05 

-48.3 

-60.8 

12.5 

4.5 

2.30 

0.5  i  0.1 

-49,0 

-59,4 

10.4 

4.0 

2.65 

0.7  t  0.1 

-51.2 

-61.7 

10.5 

3.1 

1.86 

0.6  i  0.1 

a.  Distance  between  evaporator  and  condenser, 

b.  Calcularted  by  Equation  (13) . 


In  a  number  of  runs  with  constant  evaporation  rate  and  condenser 
temperature,  the  adr  pressure  wais  increased  in  a  few  steps  from  10“3  to  0.4 
millimeters  of  mercury-  This  caused  the  evaporator  temperature  to  increase 
correspondingly.  From  temperature  and  pressure  measurements,  P  and  TTcq 


were  derived;  could  be  calculated  from  Equation  (14);  w  =  Wq  A-  TT  J 

,  TTes/ 

Only  for  high  air  pressures  was  the  difference  between  TTes  "IT  insig¬ 
nificant,  Thus,  from  each  run,  two  to  six  values  of  (P-7T')/(P-  7i^g)  were 
obtained,  ^ 
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. , .  ^  following  equation  that,  for  a  constant 

^ublimation  rate  and  average  absolute  temperature  f ,  this  £om  has  tc  be 
^d^endent  of  pressure.  So  for  each  run  the  average  value  ?P-  TT'^  / 
was  calculated  as  shown  in  Table  II.  ~  ' 


TABLE  II.  SUBLIKATION  RATES  rOR  ICE  RHERE  PRESSURE  INCREASED 
FROM  10"'^  TO  0.4  MILLIMETERS  OF  MERCURT 


l,inm 


20 


44 


w  exp-^lO^ 
gm/cm^sec  tcs, 


oc 


Vo" 


3/2  M 


RTo*"/^  Iw 


P  -  TT'  ^ 

5/ 

^t-TT’  \ 

\P  -  TTc^ 

calc 

-TTosj 

ejqpt 


0.48 

-20 

254 

0.41 

0.66 

0.68  i  0.03 

0.48 

-35 

241 

0.42 

0.66 

0.65  ±  0.02 

0.60 

-50 

229 

0.54 

0.58 

0.56  ±  0.02 

0.76 

-39 

238 

0.67 

0.51 

0.60  ±  0.08 

0.84 

-37 

243 

0.74 

0.48 

0.49  ±  0.08 

0.95 

-36 

240 

0.84 

0.43 

0.48  i  0.08 

1.06 

-48 

235 

0.94 

0.39 

0.29  ±  0.04 

1.42 

—36 

241 

1.25 

0.29 

0.49  ±  0.10 

1.68 

-44 

236 

1.49 

0.23 

0.22  i  0.06 

0.57 

-46 

233 

1.01 

0,36 

0.40  ±  0.06 

1.05 

-47 

235 

1.98 

0.14 

0.15  i  0.05 

2.05 

-44 

238 

3.60 

0.03 

C.17  i  :.C5 

a.  Calculated  by  Equation  (17), 


For  coiqjariaon  with 

«»p;„/p-7rcsi 


Rfl 


p  -  TT' 


(15a) 


use  has  been  made  of  the  approximate  proportionality  between  the  diffusion 
coefficient  D  and  T^/^/P,  thus 


_  DqPo  r=~ 
f  VT 


(16) 


Subscript  refers  to  values  at  0®C  and  one  atmosphere.  Substituting 
Equation  (16)  becomes: 

-RTq  3/2  1  w 

VoM  7^, 


P  -  //  • 


P  -  TT. 


cs/ 


=  expt 


(17) 
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In  Equation  (17)  w  is  the  mass  flow  of  water  vjq)or  per  unit 
surface  of  the  diffusion  space. 

In  the  experimental  arrangement,  the  surface  available  for  diffusion 
was  circulau*,  75  millimeters  (mm)  in  diameter,  a  distance  of  about  six  mm 
from  the  evaporator  surface.  For  the  rest  of  the  distance,  1,  a  diameter 
of  96  mm  was  available.  From  this,  the  average  flow  density  w  can  be  cal¬ 
culated  to  be  0.73  Wejrot  i'or  1  =  20  mm  and  0.66  Wgxpt  for  1  *  44  mm.  These 
values  were  used  in  the  computations  for  Table  II. 

Bradish§i£/  quotes  Knudsen  that  the  maximum  rate  of  vaporation 
frcm  an  ice  surface  at  absolute  temperature  T  is  expressed  by  the  equation: 

^max  *  Qe244(^  gm/cm^sec  (18) 

{T 

where  *  equilibrium  vapor  pressure  of  ice  in  mm  Hg  (sat.) 

0.244  *  function  of  molecular  weight  of  vaporizing  material 

(f  *  condensation  coefficient,  which  expresses  fraction  of  vapor 
molecules  entering  mass  of  ice,  following  collision  with  its 
stuface 

T  *  absolute  temperature 

Bradish  accepts  (f  equal  to  one.  Implication  of  (/equal  to  one  is  t*  i* 
every  molecule  colliding  with  the  ice  surface  is  immediately  condensed, 
although  it  may  be  evaporated  later.  Thus,  every  molecule  leaving  the 
evaporating  surfaces  will  suffer  collisions  in  the  partially  dry  material 
and  interspace  and  will  eventually  collide  with  and  condense  on  the  surface 
of  the  frozen  material  or  the  condenser  ice.  Effective  rate  of  sublimation 
is  not  equal  to  the  rate  of  vaporization  from  the  surface  of  the  frozen 
material,  bu'l  to  thr'  mass  of  water  vapor  actually  transferred  in  unit  time 
from  frozen  material  to  condenser  ice. 


C.  HEAT  TRANSFER  DJ  FREEZE  DRYING 

Lambert^/ gives  the  following  equation  for  conductive  heat  transfer  in 
freeze  drying. 


(19) 
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where  m  =  drying  rate  per  unit  area 

h  *  heat  transfer  coefficient  of  heating  mediua 
B  =  thiclmess  of  containing  wall 

X  distance  from  containing  wall  to  subliming  interface 
Ei  »  resistance  at  container  ice  boundary 
ty**  tes^erature  of  heating  medium 
kjj  *  thermal  conductivity  of  container 
tg  s  tea^erature  of  subliming  surface 
As  *  heat  of  sublimation 

k£  ■  thermal  conductivity  at  container  ice  boundary 

As  the  ice  is  co(q>letely  sublimed  the  tea^erature  of  the  product 
rises.  According  to  Abelow  and  Flosdorf,12/  this  represents  maximum 
dehydration  as  far  as  is  possible  at  this  stage  of  the  cycle.  Khen  the 
product  reaches  a  properly  chosen  control  tes^erature  the  heat  input  is 
reduced  by  alloiring  the  platen  tes^erature  to  drop.  The  product  tessera- 
ture  is  allowed  to  rise  to  a  pimdetermined  balancing  ten^jeratura  cqtud  to 
the  final  ten^jeratu.re  of  the  platen.  The  product  is  retained  at  this 
teinQjeratxure  over  the  so-called  bzdancing-out  period. 

Abelow  and  FloadorfiO/  concluded  that,  compared  with  the  flat- 
tray  design,  ribbed-tray  design  with  a  given  loading  yields  greater  out¬ 
put  at  any  usable  teoperature  condition  tested.  In  addition,  ribbed 
design,  compared  with  flat  design,  shows  more  rapid  increase  in  output 
as  heating  intensity  is  increased.  To  coiiq)ensate  for  volume  occt^ied 
by  the  ribs,  the  Icyer  of  product  should  be  thicker  to  maintain  comparable 
tray  loading  based  on  drying  platen  area.  Actual  metallic  tray  heating 
surface  to  which  pellets  were  e:qposed  in  ribbed  trays  is  ^proximately 
215  per  cent  greater  than  in  flat  trays. 

The  100®  to  40°C  heat  range  with  black  ribbed  trays  produces  the 
best  output  with  acceptable  fined,  content  of  residual  moisture. 

In  order  to  calculate  the  rate  of  drying  during  the  constant-rate 
period,  either  the  mass  transfer  equation  or  the  heat  transfer  equation 
may  be  used: 

Mass  Transfer:  w  *  k'y  (Hi-H)A 


hy  (t-ti)A 


Heat  Transfer:  w 


(20) 
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wliere  w  =  rate  of  evaporation,  pounds  per  hour 
A  =  diying  area,  square  feet 

hy  =  heat  transfer  coefficient,  Btu/Cft^)  (hr)  (op) 

k'y  *  mass  transfer  coefficient,  lb/(ft2)  (hr)  (unit  humidity 
difference) 

H£  =  humidity  of  air  at  interface,  pounds  H20/pounds  dry  air 
H  =  humidity  of  air,  pounds  H2O  per  pounds  dry  air 
t  *  ten^jerature  of  air,  °F 
t£  *  temperature  at  interface,  °F 
Ai  *  latent  heat  at  temperature,  tj[,  Btu  per  pound 

When  the  air  is  flowing  paurallel  with  the  surface  of  the  solid,  the 
heat  transfer  coefficient  may  be  estimated  by  the  dimensional  equation: 

hy  -  0.128 

where  G  is  the  mass  velocity  in  lb/(ft2)  (hr) 

When  the  flow  is  perpendicular  to  the  surface,  the  equation  is: 

.  hy  -  0.37 

D.  INTERPHASE  JUSS  TRANSFER 


From  Schrage,®/  at  the  evaporator; 
w/wg*  = 


V  =  Cfds 


Po  1 

T 

1/2 

1  - 

0 

P  * 
l^S 

^0 

\ 

\  i 

* 

*  T 

p  -  p  (_: 


1/2 


where 


MW 


Tsi 

Pounds 


PTfRTg*  (inmHg)  (ft2)  (hr) 


(21) 

(21a) 
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2447.5 

For  water,  \/Oo\  in  pounds  per  millimeter  of  mercury- 

I'  SM  }  square  foot-hour 

For  the  condenser: 


fg*  ~  ten5)erature  and  equilibrium  vapor  pressure  of  ice, 
at  condenser  interface 


Pc»  "^c  “  interfacial  vadues  in  gas  phase. 
Accurate  Approximation: 


Cf‘l 


1/2 


1 


E.  MASS  TR.4NSFER  THROUGH  DRIED  MATERIAL 

Movement  of  water  through  a  partiaUy  dry  cake  of  material  above  the 
subliming  interface  is  attributable  to  one  of  the  following; 

(a)  Bulk  flow  under  a  pressure  gradient 

(b)  Diffusion  o^  water  vapor  in  gas  phase 

(c)  Liquid  diffusion  through  solid 

During  the  constant-rate  period,  drying  usually  occurs  from  the  receding 
ice  zone;  during  the  falling-rate  period  drying  proceeds  through  the 
dried  cake. 

The  actual  structure  of  the  dried  layer  is  affected  by  solids  content, 
type  of  material  being  dried,  rate  of  freezing,  solids  redistribution 
during  freezing,  and  shrinkage  due  to  thawing  during  drying, 

Strickland-Constable  eind  Bruceii/  developed  the  following  equation 
for  the  flow  water  vapor  through  freeze-dried  clay: 
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_  (Ap)^ 

■  HnL  (23 

where  w  =  drying  rate 

L  =  thictoess  of  dry  layer 

Km  =  specific  resistance  of  the  material  being  dried 

Ap  =  pressure  drop  through  the  cake. 

As  the  downstream  pressure  is  negligible,  A  P  is  equal  to  the  equilibrium 
vapor  pressure  at  the  temperature  of  the  subliming  surface.  The  above 
relation  has  doubtful  application. 

Stephensonii/  <md  CarmanM''^  combined  classical  diffusion  and  kinetic 
theories  to  characterize  mass  transfer  in  vacuum  distillation  or  sublima¬ 
tion.  Three  operating  regions  were  identified  by  ratios  of  eqi^ibrium 
vapor  pressure  at  subliming  interface,  Pg*,  to  total  pressure,  7r» 

Region  Ps*/7T  Rate-Controlling  Step 


<^1.0 


yl.O 


^.0 


diffusion  between  evaporator  eind  condenser 

molecular  sublimation,  interphase  mass  transfer 
control 1 ing 

trauisition  region,  neither  process  controlling 


Linear  velocities  of  water  vapor  may  be  sonic  in  pores  of  dried  material 
above  subliming  surface,  if  the  resistance  to  vapor  transport  in  these 
pores  is  high.  Thus,  when  the  critical  pressure  ratio  for  water  vapor  of 
0,55  is  reached,  a  further  decrease  in  the  downstream  (space  between 
evaporator  and  condenser)  pressure  has  no  effect  on  drying  rate  or  sub¬ 
liming  surface  temperature.  This  effect  has  been  observed  by  Flosdorf.M/ 

The  Grashof  Xumbgr  used  for  correlation  of  heat  and 

mass  transfer  data:—' 


where  g  =  acceleration  of  gravity 

L  =  chairacteristic  length  of  system 
U  =  kinematic  viscosity 
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M  =  molecular  weight 
T  =  ten5)erature 
subscripts: 

Q  =  evaporator 
c  =  condenser 
h  =  heat 
m  =  mass 

Jacob  indicates  that  with  heat  transfer  Nqj.  as  low  as  1600,  free 
convectional  contribution  to  heat  transfer  in  an  enclosed  horizontid 
air  layer  is  rjgligible.  Therefore,  it  seems  unlikely  that  free  con¬ 
vection  mass  iii^sfer  is  a  significant  factor  in  freeze  drying,  as  Ngj. 
is  almost  always  below  1100. 

Guthrie  and  tfakerlin^®/  show  that  turbulent  flow  is  rarely  a  con¬ 
sideration  in  vacuum  work  below  2  mm  of  Hg,  The  same  authors  give  the 
following  criteria  which  define  the  free  molecular  and  viscous  regions 
in  terms  of  the  Knudsen  number: 

Free  molecular  flow  (turbulent)  Ni&i  ^0.33 

Transitional  region  iO  <0.33 

Laminar  flow  (viscous)  <10"2 

where  Njjjj  =  1/L 

1  =  mean  free  path 

L  =  shortest  dimension  characteristic  of  flow  cliannel 

i  ~t/p 

In  commercial  freeze  dryers  the  flow  from  the  evaporator  to  the  con¬ 
denser  is  almost  always  lamineir.  Flow  through  partially  dry  solid  above 
the  subliming  surface  is  usuailly  free  molecular  or  transitioned. 

Bulk  flow  under  a  pressure  gradient  through  a  porous  cake  may  be 
free  molecular,  transitional,  or  laminar;  although  free  path  considerations 
indicate  that  free  molecular  usually  predominates.  For  free  molecidar 
flow,  collisions  between  gas  and  solid  molecules  are  the  most  frequent, 
and  flow  of  water  vapor  is  independent  of  air  pn'ssure.  Drying  probably 
occurs  by  process  of  successive  evaporation  and  condensation. 
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AssuminE  drying  takes  place  only  at  the  free  ice  surface,  Guthrie  and 
Hakerlingi^'  present  a  relation  for  free  molecular  flovf  through  a  channel 
connecting  two  regions  of  different  pressure  derived  from  the  theory  of 
Equation  (5).  This  e3q)ression  is  adaptable  to  the  flow  of  water  vapor 
through  a  porous  dry  layer  if  a  suitable  dimensionless  constant  (correction 
factor)  Kfjj,  dependent  on  flow  geometry  and  void  fraction  of  material  being 
dried,  is  added.  It  is  further  assumed  that  partial  pressure  of  H2O  varies 
continuously  with  cake  thickness: 


'^bf  =  Kfa 


(25) 


W5£*  =  bulk  flow  rate,  mass/time-area 

d  =  average  diameter  of  flow  chemnel  (ice  crystal,  diameter) 

=  defined  in  Equation  (21a) 

For  the  transitional  flow  region,  I&iudsenil/  has  derived  the  formula 
for  pressure  drop  in  long  circular  pipes.  This  formula  is  anzd-ogous  to 
Equation  (25)  except  for  a  multiplying  factor,  J(>tei)»  Guthrie  and 
Ifzikerlingi®/  have  tabulated  this  factor  for  various  Khudsen  numbers. 

By  introducing  this  term,  Equation  (25)  may  be  applied  to  the  present 
case. 

Intennolecular  gas  phase  collisions  are  most  common  in  the  viscous  or 
laminar  region,  and  flow  rats  is  rot  independent  of  the  air  pressure.  7‘ct. 
starting  point  for  most  investigations  cf  flow  through  porous  meuia  is 
Darcy’s  Law; 


where  cake  permeatbility  in  units  of  (length)^ 

/o  =  average  gas  density  in  the  cake 
gas  viscosity 

The  right-hand  side  of  this  equation  is  approximately  correct,  since  the 
partial  pressure  of  the  air  in  the  cake  is  usuailly  small  compared  with  the 
partial,  pressure  of  the  water  v^or.  Generally  speadcing,  the  permeability 
is  proportional,  to  the  square  of  the  effective  flow  channel  diameter,  d, 
and  a  constant,  Kl,  dependent  on  flow  geometry  aind  fraction  voids. 

The  diffusion  of  water  vapor  in  the  gas  phase  occurs  only  when  Nj^  is 
in  the  laminair  flow  region; 
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Wd  =  Dv 


3  X 


(27) 


Wjj  =  rate  of  diffusional  mass  transfer 
Dy  =  diffusion  coefficient  for  water  vapor  in  air 
=  water  vapor  density  =  PM/RT 
C -  void  fraction 


Equations  (26)  and  (27)  give  the  total  mass  flux  of  water  vapor  at  any 
position  in  the  dry  cake  for  laminar  flow. 

During  constant-rate  periods,  the  partial  pressure  gradient  may  be 
assumed  independent  of  time  (steady  state).  Integration  of  the  individual 
flow  equations  between  L  and  x,  the  boundaries  of  the  dry  layer,  yields 
the  net  flow  rate  if  no  change  in  flow  mechanism  occurs  between  the  limits: 


Flow,  molecular:  w 


Po-^i 

L-x 


Viscous 


w 


Pp-Pj 

L-x 


Transitional 


^o~Pi 

L-x 


L 


0s  ^  ^  (N’Kn)  J 


(28) 

(28a) 


where  Pq  and  Pj_  are  the  partial  pressure  of  water  at  the  free  ice  and  dry 
cake  surfaces,  respectively. 

After  the  free  5ce  is  removed,  liquid  diffusion  mechanism  must  be 
considered  and  nonisothemal  unsteady  state  conditions  are  encountered 
throughout  the  porous  bed.  The  rate  of  flow  per  unit  area  by  liquid 
diffusion,  Kp,  is: 

"D  -Dl,  -|l  <1  -  f  (5  (: 

where  Dl  =  diffusion  coefficient  for  water  through  the  solid 
w  =  moisture  content  per  unit  mass  of  dry  solid  at  x. 

(5  =  bulk  density  of  dry  solid 
C  =  void  fraction 
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The  tenn  containing  the  void  fraction  accounts  for  the  fact  that  the  area 
for  liquid  diffusion  is  less  than  the  cake  cross-sectional  area.  Addition 
of  Equation  (29)  to  the  corresponding  equation  for  vapor  diffusion  or  bulk 
flow  gives  the  total  rate  of  moisture  flow  at  x.  Thus,  for  free  molecular 
bulk  flow,  total  rate  of  moisture  movement,  Wx,  is: 


w. 


ap 


-  (1  -  O' 


(30) 


A  water  balance  on  thickness  element,  dx,  of  the  slab  is  obtained  by 
equating  the  difference  between  the  input  and  the  output  of  water  to  the 
change  in  moisture  content  within  dx.  The  difference  between  the  input 
emd  output  water  is  found  by  differentiating  Equation  (30)  with  respect 
to  X  and  multiplying  by  dx. 


'^in  “  >’out  =  dx 

3x 


(30a) 


Water  content  of  the  element  is  made  up  of  the  sum  of  the  liquid  and 
vapor  portions: 


(a)  Liquid  water  in  the  solid/unit  cake  area  is  equal  to  w0  dx 

(b)  Water  vapor  in  the  pores  of  the  solid/unit  area  is  equal  to 


I  w 


c 

V. 


dx 


RIW  Cdx 
PJ 


The  change  in  the  moisture  content  of  the  element  dx  per  unit  cake  area  is 
obtained  by  di.fferentiziting  the  total  water  content  with  respect  to  time. 
The  latter  is  then  equated  to  the  input-output  of  water: 


-JL  (W6  +  Pwe)  (31) 

de 


Equation  (31)  is  a  general  equation  describing  the  mass  transfer  in 
the  falling-rate  period  of  freeze  drying.  It  may  be  simplified  if  the 
values  of  Kf^j,  /Og,  d,  Dl,  and  Q  are  assumed  constant  with  x. 


*  (1  -  €)2  6i)i, 


(32) 


Graphical  solutions  to  similar  esqpressions  have  been  obtained  for  several 
boundzuy  conditions. 


F.  MASS  TRANSFER  THROUGH  THE  VAPOR  SPACE 


Bradish®/  states  that  the  retardation  effect  of  the  porous  dry  aass 
above  the  receding  ice  surface  of  the  material  being  dried  is  small.  He 
believes  that  the  transfer  of  vapor  is  limited  by  the  rate  of  diffusion. 

Bradish§/  theorizes  that  the  interspace  gap  is  so  short  that  molecules 
are  transferred  almost  without  collisionj  then  the  sublimation  rate  is 
equal  to  the  rate  of  vaporization  at  the  surface  of  frozen  material  minus 
the  rate  oi  vaporization  at  the  surface  of  the  condenser  ice. 

If  ntmerous  collisions  occur  in  the  interspace,  this  siaq)le  difference 
equation  no  longer  applies.  It  must  include  terms  expressing  proportion  of 
molecules  which  traverse  the  interspace,  as  distinct  from  those  which 
eventually  condense  on  the  surface  from  which  they  were  evaporated. 

In  drying  systems  in  idiich  evaporation  and  condensing  surfaces  are 
widely  separated,  or  in  which  the  permanent  gas  partial  pressure  is  high, 
the  subnotion  rate  is  low  as  a  consequence  of  the  numerous  molecular 
collisions  in  the  interspace. 

Assuming  parallel  condensing  and  evaporating  surfaces,  the  flow  channel 
between  the  sa]iq)le  and  the  condenser  is  a  rectangular  duct.  If  free 
molecular  flow  pertains,  the  rate  of  water  vapor  transport,  w,  is: 


where  a,b  ■  dimensions  of  rectangular  cross  section 
Z  «  distance  from  condenser  to  cake 

K  *  constant,  dependent  on  ratio  a/b  (tabulated  by  Guthrie 
and  Kakerling  when  a/b  is  l/2,  K  is  1.151) 

In  the  viscous  region  (Nrh  less  than  0.01),  the  va?)or  movement  is  controlled 
either  by  vapor  diffusion  through  stagnant  air  or  by  diffusion  and  bulk 
flow  under  a  pressure  gradient.  Stefan’s  diffusion  equation  described 
the  former  case: 

fHW  (Pj-Pc) 

KlZpin  (34) 

Dv  =  diffusion  coefficient  for  water  vapor  through  air 
HR  =  molecular  weight  of  water 
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Pljji  =  log  mean  air  pressure  =  pg-p^/ln  Pg/Pi 
IT  =  total  pressure  (p^  +  p^) 

Z  *  distance  from  cake  surface  to  condenser 


P  refers  to  vater  vapor,  p  to  air 

For  ideal  gases  at  constant  T  and  total  P,  a  general  diffusion  equation 
(neglecting  themed,  pressure  and  forced  diffusion)  is: 

if  Ml  Dr  MK  dP 

7  *  M  Ma  ’  ST  d^  (35) 

V  *  mass  flux  of  water  vapor  per  unit  time-  and  area 
Wa  “  mass  flux  of  air  per  unit  time  and  eurea 
P  *  partial  pressure  of  water  vapor  at 
^  »  distance  from  sample  surface  to  condensing  interface 
Ma  *  molecular  weight  of  air 

An  air  material  balance  around  the  dryer  shown  below  allows  estimation  of 
Wa  in  terns  of  speed  of  vacuum  pvmp,  S,  and  volume  rate,  L,  at  which  air  is 
leaked  into  the  space  bounded  by  dotted  line: 


RT 


oo 


*J%^ 


A3 sumptions: 

(a)  Steady  state  operation  (w  is  equal  to  K). 

(b)  Gas  composition  at  any  cross  section  is  constant. 

(c)  Air  only  is  removed  by  vacuum  pump  (condenser  and  trap  remove  all 
water  vapor) . 

(d)  Total  pressure  of  system  is  insignificant.  _ 

Therefore,  air  flow  at  ^  is: 

TrMa 

Wa  =  (S  -  L)  art  where  A  =  area,  (36) 

Substituting  this  in  Equation  (35) 


i-  ,  MW  (S  -  L)  (TT  -  P)  ^  Dy  MW  ^ 

IT  ^  "  ART  ^  RT  ^  (36a) 


and 


dp  dp 

d^  d 


stag,  air 


,  (TT  -  P)  MW  (S  ■  L) 
wARTp 


(36b) 


Equation  (36b)  b«  integrated  from  p^  (  ^  *  0)  to  p^  (  ^  =  Z)  to  obtain 

(37) 


w 


Dv  MWTr^ 

1 

(rT.  j 

1  Pin 

9  1 

MW  (S  ■  L)  7r 
ART 


In  this  equation 


Pc-Pi 


liRTAwpc  -  MW  (S  -  L)  TIP 


[RTAwpi  -  MW  (S  -  L)  TT 


'glc 
IT  Pi 


or  a  modified  log  mean 
air  pressure.  (38) 


When  (S  -  L)  =0,  Equation  (37)  reduces  to  the  equation  for  diffusion 
through  stagnEmt  air.  Equation  (34). 


G.  OVER-ALL  MASS  TRANSFER  EQUATION 

Pp*  Pjf_  (S  -  D-  Zpja 

r  s  ~  Pc  ~ 

w  =  “l  T^x  .  RTZPiffl 

^  TIA  5OT  55^ 


(39) 
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Assumptions: 

(a)  Equations  (21)  and  (21a)  apply  to  mass  transfer  at  subliming  and 
condensing  interface,  respectively 

(b)  Equation  (28)  applies  to  transport  throu^  dried  cake 

(c)  Equation  (37)  applies  to  vapor  flow  from  dry  cake  to  condenser 

General  Form:  w  =  driving  force  /  sum  of  resistance 
Denominator  includes  resistances  at  two  phase  interfaces,  through  dried 
cake  and  between  evaporator  and  condenser.  Free  molecular  flow  througji 
cake  layer  and  laminar  flow  to  condenser  have  been  assumed. 


H.  OVER-ALL  MSS  TRANSFER 

Table  III  gives  a  svunmary  of  mass  transfer  resistances  for  all  possible 
mechanisms  of  vapor  movement.  The  relative  importance  of  interphase  mass 
transfer  resistances  in  pure  ice  sublimation  is  shown  for  various  totsd. 
pressures  and  drying  rates  in  Figure  5.  This  plot  is  based  on  theoretical 
calculations  using  Equations  (21),  (22),  and  (34)  and  assuming  equilibrium 
water  v8q)or  pressure  at  condenser  of  30  microns  (Tg*  =  -50°C). 


TABLE  III,  MASS  TRANSFER  RESISTANCES  FOR  POSSIBLE 
MECHANISMS  OF  VAPOR  MOVEMENTS' 


RESISTANCES 

FREE  MOLECULAR 

MASS  TRANSFER  STEP 

OR  TURBULENT 

LAMINAR  OR  VISCOUS 

Interphase  I 

1 

1 

(ice  to  water  vapor) 

\^s  0s 

Ps/^3 

Dry  cake 

(Subliming  surface  to 
cake  surface) 

L-x 

Kfn  /%d 

L-X 

DvMWe^,, 

RT  M 

Evaporator  to 
condenser  (cake 
surface  to  condenser) 

(a+b)Z 

RTZplm 

8/3  Kab^ 

HK  D^TT 

Interphase  II 

X 

1 

(water  vapor  to  ice) 

Transfer  Rat.  As  Per  Cent  of  Total  R-^t. 
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Figure  5. 
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u  Tr^msfer  Resistance  in  Sublimation  of  Ice  as  a 
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Combination  of  Equation  (19)  for  heat  transfer  and  Equation  (39)  for 
mass  transfer  aiUows  specification  of  temperatures  in  the  ice  cake  during 
the  constant-rate  period,  since  tg  (or  T_*)  and  Pg‘“‘  are-em- equilibrium 
pair.  Effect  of  dry  cake  mass  transfer  resistance  on  subliming  inter¬ 
face  temperature,  is  shown  in  Figure  6. 

Ttie  possible  freeze-drying  mass  transfer  driving  forces  and  resistances 

arej 


Driving  Force 


Resistance 


?fat2r  vapor  partial  pressure  dif¬ 
ference  between  evaporator  and 
condenser 

Bulk  flow  of  gas  to  pumping 
Transport  by  natural  convection 


Entrainment  of  ice  particles 
by  moving  gas  stream 


Resistance  to  mass  transport  at 
the  subliming  interphase 


Diffusional  or  flow  resistance 
through  dry  cake  above  interface. 

Diffusional  resistance  to  water 
vapor  flow  afforded  by  inert  gas 
between  evaporator  and  condenser 

Resistance  to  interphase  transfer 
at  the  condenser 


Thermal  diffusion 


Subliming  Surface  Temperature, 


Figure  6.  Diy-Cake  Mass  Trzmsfer  Resistance  as  a  Function  of  Subliming 
Interface  Temperature  at  Various  Total  Pressures  and  Hass 
Transfer  Rates. 
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III.  SAMPLE  CALCULATION 


Given  a  freeze  dryer  heated  by  conduction.  The  condenser-to-evap orator 
spacing  is  0.5  foot.  The  evaporator  area  is  20  square  feet.  The  refrig¬ 
erated  plate  condensers  which  are  maintained  at  -40°F  are  used  to  dry  a 
one-.jich-thick  slab.  The  cake  is  in  perfect  thenaal  contact  with  the  tray 
bottom..  The  material  thaws  at  +2°F.  The  maximum  allowable  tenq)erature  is 
0°F. 

Calculate: 

1,  Maximum  possible  drying  rate  and  heat  input  if  P  =  1000 

2.  Csdctilate  tg  or  Tg->  and  ~fl' if  the  total  heat  input  is  set  at  16,300 
Btu’s  per  hour. 

Neglect  the  dry  cake  resistance  and  assume  that  the  thermal  conductivity  of 

Btu 

the  material  is  1.5  •  The  relationship  between  vapor  pressure  and 

teD5)erature  is  the  same  as  that  for  water,  and  the  diyer  efficiency  is  75 
per  cent.  Assume  that  the  ice  layer  recedes  uniformly  into  the  cake. 


Solution  to  1, 


Assume 


As  H  “  150  Btu/hr  ft^ 


As  11 

t.-t  . -  =150x-^x  ^ 

^  s  Ki  12 


1.5 


8.3°F 


tg  =  -  8.3°r 

The  corresponding  drying  rate  is  0,130  Ib/hr  ft^  or  the  calculated  value 
of  AsK  =  (  A  s  calc  =  0.130  x  1220  =  156  Btu/hr  ft^ .  Since  the  first 
approximation  was  150  Btu/hr  ft^,  we. get 

w  =  0.130  Ib/hr  ft^  =  2.5  Ib/hr. 

The  true  amount  of  heat  added  is 

156  Btu/lir  ft^  =  208  Btu  -per  hr  ft2 
0.75 


or 


208  Btu 


hr  ft^ 


=  X  20  ft2  =  4160 


Btu/hr , 


Solution  to  2. 


The  value  of  the  initial  evaporating  surface  ten5)erature,  ts,  is  agadn 
found  using 


Ki 


+  t. 


ts 


0.75  X  16^300 
1.5  X  12  X  20 


+  0  =  -340F 


The  vapor  pressure  of  water  at  the  subliming  interface,  Pq,  is  found 


Po  = 


P  * 


1  - 


Vs’ 


lb 


^3  "  mm  Hg  ft^  hr 

*  0.140  mm  Hg 

Kg*  «  Pg*  16.61  lb/ft2  hr 

1220  x' 26  ’  Ib/ft^  hr 


s  r  s 

Eq  0.75  X  16,300 


Eg 

fg  «  0.982 


Pp  =  0.140  X  (l-0.030)/0.982  =  0.138  mm  Hg. 


Similarly,  the  partial  pressure  of  water  at  the  condenser  interface 
is  computed  to  be  0.104  mm  Hg.  The  total  pressure  is  then  calculated, 
neglecting  bulk  flow. 


Rearranging, 

426  2.334 

(0.8525)  (^)  (18) (0.138-0.104) 


Im 


(0.7302)(426)(0.5)(0.5) 


=  0.005  mm  Hg 


where  the  value  of  T  has  been  taken  as  equal  to  the  evaporator  interface 
teii?)erature .  The  totail  pressure  is  calculated  from  the  log  mear  air 
pressure  as  follows: 
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Pq-Pq 
In  7r-  Pc 
TT-  Po 


0.138  -  0.104 

In  -  0»104 
TT-  0.138 


=  0.005  mm  ilg 


When  this  equation  is  solved  for  IT  ,  we  get  IT  =  -0.138  mm  Hg. 

Calculation  of  Drying  Time 

Calculate  the  drying  time  for  the  material  in  the  above  problem,  if  the 
criticail  moisture  is  known  to  be  1.5  lb  H20/lb  solid.  Drying  tests  indicate 
that  the  rate  is  proportional  to  the  moisture  content  for  the  entire  falling- 
rate  period.  The  initial  solids  content  is  10  per  cent;  the  bulk  density  of 
product  is  six  Ib/ft^,  and  the  final  moisture  of  the  product  is  one  per  cent. 

^  ^  Eq  a  Kf 

where  ’’a"  represents  the  slope  of  the  rate  ^  moisture  curve 
Slope  "a"  is  estimated  by  the  relation 

Wc  =  ^  =  a  ifc 

As 


or  Eq  0.75  x  16,300  0.33  lb  solid 

- - 1 -  _  - jj - 

As  Kc  1220  X  1.5  X  2.0  ft^  hr 


After  substitution  into  the  equation  for  the  total  drying  time,  v,-e  obtain: 

(1)(6)  (1220)(9.0-1.5)(20) 

6*  =  (12)  (0.75)(16,300) 

These  sample  calculations  illustrate  the  use  of  the  equations  developed 
earlier,  particulaily  Equations  (19)  and  (39).  The  heat  and  mass  transfer 
for  conductive  freeze-drying  and  the  resistances  connected  with  this  type 
are  considered  in  these  relationships.  The  theoretical  rates  and  times 
determined  here  were  subsequently  verified  e:qperimentally  and  good  agree¬ 
ment  was  obtained. 
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IV.  CONCLUSIONS 


The  principles  of  freeze  drying  are  relatively  recent  in  origin  and 
are  published  in  a  variety  of  research  areas.  This  report  presents  those 
principles  governing  freeze-diying  rates  of  heat  and  mass  transfer,  and 
defines  them  in  practical  terms. 

Some  of  these  principles  are: 

(a)  Primary  factors  which  determine  subliming  interface  temperatures 
are  absolute  total  pressure  and  the  mass  transfer  resistance  of  the  dry 
solid  above  the  evaporating  surface.  During  the  drying  of  substances 
containing  solids  the  equilibrivm  vapor  pressure  is  usually  greater  than 
the  total  pressure,  but  approaches  it  as  the  resistance  to  vapor  flow 
through  the  residual  solid  above  the  interface  decreases. 

(b)  The  rate  of  drying  is  determined  by  the  dryer  efficiency  and  heat 
input  to  the  evaporator.  ?Ihen  drying  slabs,  the  heat  of  sublimation  must 
be  conducted  through  the  frozen  layer.  Therefore,  heat  input  and  slab 
thickness  govern  the  rate  of  drying;  the  total  pressure  and  conductivity 
determine  the  drying  temperature  level. 

(c)  The  drying  time  per  unit  mass  is  greater  for  pellets  than  for 
slabs.  The  temperature  control  in  a  bed  of  pellets  is  more  difficult. 
Even  though  drying  occurs  throu^out  the  pellet  bed,  these  two  factors 
are  related.  Since  the  outer  surface  of  any  pellet  dries  first,  the  heat 
necessary  for  sublimation  during  later  stages  must  pass  through  this  dry 
exterior.  Thermocouples  tend  to  yield  erroneous  measurements  tmless 
properly  located,  and  the  heat  input  must  be  reduced  earlier  to  avoid 
overdrying  of  the  already  dry  exterior  of  any  pellet,  thus  increasing 
the  drying  time.  These  same  considerations  preclude  the  use  of  radiant 
heat  in  a  freeze-drying  operation,  because  penetration  of  radiant  energy 
into  frozen  sampleis  is  small  and  heat  inputs  must  be  kept  low  to  prevent 
overdrying  or  uneven  drying  of  the  materials. 
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GLOSSARY 


All  basic  symbols,  subscripts  and  superscripts  used  in  this  memorandum 
are  defined  here  except  for  a  few  terms  used  only  in  limited  contexts. 

The  latter  are  defined  in  conjunction  with  thej.r  use.  Where  a  symbol  is 
shown  as  having  more  than  one  meaning,  the  context  will  indicate  which 
definition  is  pertinent.  The  primary  dimensions  listed  have  the  following 
meanings: 


L  =  length 

M  =  mass 

0  =  time 

T  =  temperature 

F  »  force 

Q  =  energy 

PRIMARY 

SYllBOL 

MEANING 

DIMENSIONS 

a 

"slope”  constant  from  rate  ^  moisture  cui*ve 

M/l20 

A 

area 

l2 

B 

wall  or  film  thickness 

L 

c 

heat  capacity 

qMAt) 

concentration  per  unit  mass 

(tmit)/^ 

G 

conduct  Iviiij 

mlA0 

D 

diffusion  coefficient 

E 

dryer  efficiency 

- 

F 

volume  flow  rate 

\?IQ 

g 

acceleration  of  gravity 

CM 

h 

heat  transfer  coefficient 

Q/L^  6^(At) 

J 

"transitional"  flow  correction 

- 

k 

thermal  conductivity 

Q/L  0(At) 

K 

"geometric"  factor  for  noncircular  flow  channel 

smoL 

MEANING 

PEIMARY 

DIHENSIONS 

1 

molecular  free  path 

L 

L 

cake  thickness 

L 

characteristic  length 

L 

Mor  MK 

molecular  weight 

M/mole 

N 

dimensionless  number 

P 

partial  pressure  of  air 

Y/\? 

P 

partial  pressure  of  water 

F/l2 

q 

rate  of  energy'  transfer  per  unit  area 

q/l^ 

Q 

heat  content  per  unit  volume 

Q/L^0 

r 

radius 

L 

outside  radius 

L 

gas  constant 

FL/T{mole) 

heat  or  mass  transfer  resistances 

l2(At)  0/Q  or 

S 

pumping  speed 

t 

condensed  phase  temperature- 

T 

T 

gas  temperature  (degrees  absolute) 

T 

u 

velocity 

L/0 

V 

volxaae 

l3 

w 

rate  of  mass  transfer  per  unit  area 

M/l20 

K 

moisture  content  per  mass  of  dry  solid 

X 

voidable  distance  in  ice  or  solid 

L 

Z 

distance  from  cake  surface  to  condenser 

L 

SYMBOL 

A 

r 

A 

u 

M 

Q 

TT 

<P 

6 

f' 

I 

a 

e 

f 

w 

t 


PfilMARY 

MEANING 

DIMENSIONS 

fvinction  defined  by  Equation  (21a) 

dfh 

difference 

function  defined  by  Equation  (9) 

latent  heat  per  unit  mass 

kM 

kinematic  viscosity 

VLfL^Q 

absolute  viscosity 

M/L0 

ratio  of  liquid  to  ice  density 

toxal  pressure 

f/l^ 

function  defined  by  Equation  (8) 

bulk  density  of  solid 

m/l^ 

true  density 

vl/\? 

accommodation  coefficient 

time 

9 

void  fraction 

cake  permeability 

variable  surface  from  cake  surface  to  gas  phase 

L 

SUBSCRIPT 


MEANING 


a 

b 

an  r* 

container 

c 

condenser 

d 

conduction 
constant  rate 
critical 
diffusion 

e 

effective 

E 

evaporator 

produced  or  generated 

f 

fusion 

h 

falling  rate 
final. 

heat  transfer 

i 

ice 

L 

ice-container  interface  gas  at  the  cake 
liquid  or  liquid  phase 

m 

mass  transfer 

0 

maximum 

dry  cake  or  material 
initial  or  datum 

R 

gas  at  the  subliming  interface  orifice 
refrigerant 

s 

solid 

V 

sensible 

sublimation  or  subliming  interface 
vapor 

X 

value  at  position  x 

w 

water 

Gr 

latent  heat 
heating  medium 

Grashof 

Kq 

Knudsen 

Nu 

Nusselt 

Pr 

Prandtl 

Re 

Reynolds 

In 

log  mean 

bf 

bulk  flow 

fm 

free  molecular 

SUPERSCRIPT 

(x) 

■tt 

FUNCTIONAL  I 


erfcy  = 


MEANING 

average  or  mean  (x) 
equilibritan  value 


-Ei(-y) 


